The probability for a wave packet to remain in a disordered cavity 
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We show that the probability that a wave packet will remain in a disordered cavity until the time t 
decreases exponentially for times shorter than the Heisenberg time and log-normally for times much 
longer than the Heisenberg time. Our result is equivalent to the known result for time-dependent 
conductance; in particular, it is independent of the dimensionality of the cavity. We perform non- 
perturbative ensemble averaging over disorder by making use of field theory. We make use of a 
one-mode approximation which also gives an interpolation formula (arccosh- normal distribution) 
for the probability to remain. We have checked that the optimal fluctuation method gives the same 
result for the particular geometry which we have chosen. We also show that the probability to re- 
main does not relate simply to the form-factor of the delay time. Finally, we give an interpretation 
of the result in terms of path integrals. 

PACS numbers: 73.23.-b, 03.65.Nk 



The interest of experimentalists in open quantum dotsS 
motivates the computation of the probability to remain 
in a weakly disordered cavity. This problem is similar to 
the computation of the time, dependent conductance of 
weakly disordered samplesoQ. However the geometry of 
our problem allows a simple computation scheme, which 
is similar to the zero-dimensional non-linear er-model. 

Let us assume that we have a weakly disordered open 
cavity and the whole system is filled by non-interacting 
fermions at zero temperature. At time t = the Fermi 
level outside of the cavity decreases and particles begin 
to escape from the cavity. The probability for the par- 
ticles to remain in the cavity of volume 5* till time t is 
given by the ratio 



Pit) 



J s g{r,t)df 
Js Q{r, 0)dr 



(1) 



where p{f, t) is the density of the particles on the Fermi 
surface, and p(r, 0) = 1 for r £ S and elsewhere. On a 
short time scale-this is the standard problem of the dif- 
fusion emissions The time evolution of the distribution 
function can be found by solving the diffusion equation 
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dt 



- DV Z 



g(?,t) = 



(2) 



with the boundary conditions nV g = at the walls of the 
cavity (ft is normal to the boundary^ and g + 0.71£nV g = 
at the open edge of the cavity.El Here £ is the mean 
free path and the scattering is uniform. In further cal- 
culations we will use the open edge boundary condition 
q = 0, assuming that £ is much smaller than all the char- 
acteristic lengths of the system. 

The solution of the diffusion equation can be repre- 
sented as a sum over diffusion modes, each of them de- 
caying exponentially at its own rate. The "lowest" diffu- 
sion mode is computed in Appendix |A| for two examples 



of circular and spherical cavities with the contact in the 
middle, see Fig. [l]. Therefore the probability to remain 
in the cavity behaves like p(t) — e (t , where 7 is the 
decay rate of the "lowest" diffusion mode. This escape 
rate is proportional to the size of the contact and may 
vary from zero to the inverse diffusion time through the 
system E c /h. In the rest of the paper we will use units 
where h = 1, and then we have in general 



< 7 < E c 



(3) 




V r g = 



FIG. 1. Boundary conditions for the diffusion mode in the 
circular cavity with the contact in the middle. 

In quantum mechanics the evolution of the density ma- 
trix is given by the product of two exact Green functions 

p(r 1 ,r 4 ,t)= J dr2dr 3 G R (f 1 ,r2,t)G A (f 3 ,f 47 t)p(f 2 ,f 3 ,0) 

(4) 
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In order to compute the probability to remain in the cav- 
ity we should prepare our system in the state 



P(r2,r 3l 0) = / dpe 



p(r 2 -r 3 ) / r 2 + r 3 



0(^-^,0), 



S(E F - H(p, 



(5) 



equation for the cavity has solutions, which are expand- 
ing waves far away from the system. The corresponding 
eigenvalues of the energy are complex E n — iT n /2, see 
Rcf. ||, §134. The delay time is the analog of the density 
of states T(E) = Im^jB - E n + iT n /2) jpd_it can 
be expressed in terms of the scattering matrixE3EJ. The 
disorder averaging leads to 



where H.(p,r) is the Hamiltonian of the system, and 
g(f, 0) is one inside the cavity and zero outside. The 
quantum-mechanical expression for the probability to re- 
main is therefore 



P(t) = 



1 



J s p(r,f,t)df _ 
J s p(r, r, 0)dr 2-kvS J 2tt 



x ( drldf 2 Gl F+f (fi,f 2 )G^_ f (f 2 ,fi) , (6) 
j s 

where v is the density of states on the Fermi surface. 
The exact Green functions depend on the realizations of 
the disorder potential. The non-perturbative averaging 
over the (5-correlated disorder potential can be done by 
making use of field theoryO. Our main result is 



p(t) = exp 



-^arccosh 2 
2A 



tA 

IT 



(7) 



where A is the mean level spacing of the "closed" cavity 
1/A = vS and 7 is the classical escape rate. This for- 
mula gives exponential decay on the short time scale and 
log - normal decay on the long time scale. Equation (0) 
was derived under the conditions 



7 < E c 
t > e: 



log(tA) < 



Lilog^a/ir)^ 
Lift 



2 dimensions, 

3 dimensions, 



(8a) 
(8b) 

(8c) 



led 



for all relation of 7 and A. Particularly, the exact prob- 
ability to remain should have some features on the time 
scale E~ Y near t — and t — 2tt/A, as it takes place 
for spectral form-factors of closed systems! Our result 
Eq. (Q) does not have any structure on the time scale 
E~ x and this is the meaning of the conditions Eqs. (8a|) 
and (|8b|). If the last inequality En. (Bq) is not fulfi. 
one should use the ballistic actionH in the field theory. 

It is interesting to compare the result Eq. (Q) with the 
form-factor of the delay time 

1 f du> u> , _ , 11), 

= 2twS J 2tt e ^(^f + -^)T(E F - — ) 

1 f duj .., , t 



2nvS / 2tt 



x /dr 1 dr 2 G| F + f (r 1 ,r 1 )G^_^(r 2 ,r 2 ), (9) 
J s 

where the subscript 7 means that our system is open 
(it has the classical escape rate 7). The Schroedinger 



Kr\t)p(t) 
tA 



^(<)=min(— ,1), 



(10) 
(11) 



where the superscript "unit" means that the form-factor 
was computed as if the system is not symmetrical under 
the time reversal, and p(t) is given by Eq. (Q). From 
the definition Eq. (^) one can see that the form factor of 
the delay time K^(t) approaches the form factor of the 
density of states K (t) when 7 goes to zero and this is 
consistent with Eqs. (0) and (|l^). The decay rate of the 
lowest diffusion mode, 7 — > 0, when the ope ning o f th e 
cavity becomes smaller L\ — > 0, see Eqs. (A3) and (|A7|). 
In this case all particles remains in the cavity forever and 
our solution gives pit) — * 1. 

Our results,— are inconsistent with the random ma- 
trix theoryc3~t5l, which predicts a power law decay of 
both Kj(t) and p{t) if the number of open channels is 
much smaller than the dimensionality of the Hamilto- 
nian. However, on the short time scale, tA <C 1, one has 
from Eq. (fy 



log(p(t)) = 7 t l 



tA 
Gtt 



(12) 



and this is similar to the numeric results of Ref. 16 and 
the random matrix theory calculations of Rcf. [17|. The 
relation between K^(t) and p(t) siiajlar to Eq. ([L0|) ap- 
pears in the random matrix theorylij too. 

Before going on to explain the averaging procedure let 
us give the semiclassical path - integral interpretation of 
the above result. Expansion of the Green function in a 
sum over classical paths j from the point f\ to F 2 is 



GE A (n,r 2 ) 



A j exp{±iS j (fi,f 2 ,E)} (13) 



where Sj (r\ , r 2 , E) is-the action along the path and Aj 
are some coefficients.113 Therefore, the probability to re- 
main is the double sum over trajectories, which can be 
separated into the diagonal and off-diagonal parts. Sta- 
tionary phase integration over uj in Eq. (0) gives 



Pdiag(i) 



A 2 

^ vS 



Sit-Tjifi^Ep)], (14) 



where Tj is the time which takes the particle with energy 
Ep to go from the point 77 to the point r 2 along the 
path j . Equation (fl4|) is equivalent to Eq. (|lj) and there- 
fore Priiag(t) = e~ 7 Ts the classical probability to remain. 
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Therefore, the log-normal tail of the quantum probabil- 
ity to remain represented in Eq. (0) is determined by the 
off-diagonal part of the sum over trajectories in Eqs. @ 
and (|13|). 

It is known that the integrals over coordinates in 
Eq. ( [l4|) can be, pomputed by using the skeleton of the 
periodic orbitaiScJ 



PdiagO) = T i 



S(t-rTj) 
^ |det(/-MJ)| 



(15) 



where the index j runs over all primitive periodic orbits 
inside the cavity in configuration space and on the Fermi 
surface in phase space. Each orbit has the period Tj and 
the monodromy matrix Mj . From Eq. ( |l5|) one concludes 
that the decay rate of Pdiag(^) is equal to the difference 
between the Lyapunov exponent and the Kolmogorov - 
Sinai entropy and should hp equal to the decay rate of 
the lowest diffusion mode. El 

The semiclassical form of the Green functions Eq. flic 
can be used for computation of the form-factor Eq, 
The stationary phase integration over r\,r%,u givea 2 ^ 



tA 
to* 



(16) 



which is in agreement with our result Eq. ( fTcf ) taken for 
short times. 

Let us assume that the leakage is so small, that all 
gradients are small on the scale of the mean free path. 
The averaging of the product of two Green functions from 
Eq. is given by the functional integral 

(G R G A ) = 2(nv) 2 [ Q^(n)k 00 Qf a (r 2 )e- F ^DQ 

JQ 2 =1 

(17) 



F[Q] = ^- S txj 



D(VQ) 2 + 2i(uj + iS)AQ 



dr 



£>rstr(VQ) 2 < 1 



(18) 
(19) 



where D is the diffusion coefficient on the Fermi surface, 
see Ref. ||, page 60, Q = A at the contact, and riVQ = 



at the walls. I 

One can use the polar decomposition of Q for calcula- 
tion of the integral. The probability to return should be 
independent of the symmetry of the system and we will 
use the simplest "unitary" case. Then the action Eq. JlS| ) 
has to be multiplied by a factor of 2 and 



Q 12 = iusin9v 1 Q 21 = -ivsinOu 



9 
%Q X 



1 - 2j]Tf 

-2r)* 



2r) 



2 VV * 

1 - 27777 
2rj* 





e ix 



~2r] 
27777* 



1 + 2kk* 
-2in* 

1 + 2kk* 
2ik* 



2lK 

1 - 2kk* 

-2iK 
1 — 2kk* 



where 9, 0±, <fi, x are commuting variables and k, k* , 77, 
77* are anticommuting variables which parameterize the 
matrix Q. Let us write explicitly the elements of the 
matrix Q which appear in the integrand in Eq. (|l7j ) 

[1 - 27777*][1 + 2KK*]e^ sin6- Ar] K *e ix sinhfli 
-2t7*[1 + 2KK*y rt > sin 9 - 2[1 + 2r)rf\K*e ix sinhfli 

— 2i[l - 2ryq*]Ke ic t> sin6» + 2it][l - 2nK*]e ix sinh6»i 
4^77* ne** sin 9 + i [1 + 27777*] [1 - 2kk*} e ix sinh 6»i 

Q 21 = -i 

( [1 + 2kk*} [1 - 27777*]e~^ sin 9 - 4 K r]*e- lx sinh 9 X 
X ^ -2m*[l - 27777*]e~^sin6»-|-2i[l-2KK*]77*e- lx sinh6'i 

-2[1 + 2KK*\r]e-^ sin 9 - 2/t[l + 27777*]e^ x sinh6»i 
+4iK*77e"'^ sin 6 + i[l - 2kk*][1 + 2r 1 'q*]e- ix sinhfli 

Our purpose is to reduce the functional integral over 
the super-matrix Q to the conventional integral. The 
minimum of the action is. reached when u and v are in- 
dependent of coordinates^! and the action becomes 



+ 2(iu> - 0)(cos6» - cosh6»i) >df, 



(20) 



where 9 = at the contact and nV9 = at the walls. 
The same boundary conditions were applied to Eq. (||) 
and therefore we can use the diffusion modes for comput- 
ing the functional integral over 9(r). Due to the condition 
Eq. (|8a|) only the lowest diffusion mode contributes and 
the functional integral becomes a conventional integral 
over the amplitude of this mode 8. The lowest diffusion 
mode is almost uniform, see Appendix [A|, and therefore 



2 df=Q 2 



i J s 9 2 df=e 2 , 

i J s D(\79) 2 dr = 1 e 2 , 



— / cosh 9dr = cosh 6 
S Js 



sinh 9dr = sinh O , 



and so on for other functions of 9 and 9\. Particularly 
the pre-exponential factor in the expression for the return 
probability becomes 
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Q\\Q{\ = -8r]ri* kk* (sm z + sinh 2 &i) 



(21) 



We computed this factor explicitly from the expressions 
for Q 12 and Q 21 . This expression was computed in Ap- 
pendix 3 of Ref. ^ but the numeric coefficients are differ- 
ent. 

The model reduces to the conventional integral 



p(t) = -2{ttu) 



! f ^_ e -iut f 8w * KK *( s i n 2 e + sinh 2 0i ) 
J 27T J y / 



_ r) 12 r ,21 



2Q21 



2 8 (cosh0! — cos©) 2 dr]drf dn* dn 



=dR ± 



/T 
— — e~ 



= J 2 =dfl 2 

_ F cosh 0i + cos© . 



cosh 0i — cos0 



sin sinh 1 d<d dQ i 



(22) 



F 



7 1— i + (iw-0) (cos 0- cosh 0! U , (23) 



where we have kept Efetov's notations for differentials 
and Jacobians, see Ref. |, page 107. 

The standard change of variables A = cos and Ai = 
cosh 0i leads to the relatively simple expression for the 
probability to remain 

m =\f * tA/ :: 2X 

2 Jmax(-l,l-tA/ir) tA/7T 



X P~ C arccos2 (A)+arccosh 2 (£A /tt+A)) 



(24) 



For both short times t <C tt/A and long times t 3> tt/A 
we can put A = 1 in the exponent on the right hand side 
of Eq. ( 24]) . Then the integral can be computed exactly 
and we arrive at our interpolation formula Eq. (j?j). The 
computation of the form-factor is similar, 



7^unit 



(*) 



1 

2 

x e 



dX 

max(— 1,1— tA/ir) 
— (arccos 2 (A)+arccosh 2 ( t A/ 7T+A)) 



(25) 



and Eq. ( |To| ) matches both short and long time asymp- 
totes of Eq. (p5|). Both results Eqs. (^4j) and ( |25| ) where 
derived under the condition Eq. (|19|) which becomes 
Eq. @ in the one-mode approximation. 

The functional integral over 0(f) in the theory with 
action Eq. (p^p^pap be computed by the optimal fluctu- 
ation methom'oH. The calculations are straightforward 
and give the same results, see Appendix |b], where we have 
checked the three dimensional case. In the same way one 
can also check the two dimensional result. 

In the present work we have assumed that the opening 
or contact is small. If it is not the case for some system 
the result for the probability to remain may be different. 
In such a system one should consider the optimal fluctu- 
ations of the random potentialc-3. In our geometry this 



method can be used for very long times, which broke the 
condition Eq. ([k]). 

In summary, we have obtained the arccosh-normal dis- 
tribution for the probability to remain in a disordered 
cavity. We made use of the one-mode approximation, 
which is similar to both the zero-dimensional er-model 
and the optimal fluctuation method. Our result is not 
consistent with the random matrix theory prediction. 
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FIG. 2. Radial dependence of the diffusion mode for 
< Li < L 2 . 



APPENDIX A: COMPUTATION OF THE 
DIFFUSION MODE. 



Let us compute the diffusion mode for the circular cav- 
ity with the circular contact in the middle, as shown in 
Fig. [I]. The lowest mode is independent of the polar an- 
gle and the radial dependence is expressed in terms of 
the Besscl functions 



e(r)ocIm H£\rM)H?\L 2 M 



= lm 



^ V D' 1 ' V D' 




(Al) 
(A2) 



where Z^y^ has to be less that the first root of Ji(x) 
and this is possible for roughly L1/L2 < 1/2. The solu- 
tion for Li <C £2 is shown in Fig. y and in this case 



D 



(A3) 
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We see that the mode is non-uniform for roughly L\ < 
r < L2/2. Therefore the presence of the contact strongly 
affects the diffusion mode in a quarter of the cavity area. 

In the case of the three dimensional sphere it might be 
difficult to maintain the contact in the middle, but it is 
still interesting to compute the lowest diffusion mode. In 
three dimensions the lowest mode is independent of polar 
and azimuthal angles and depends only on the distance 
r from the center. The lowest mode in three dimensions 



, 1 • lor (r-Li) 2 
g(r) oc - sm a SLi j-^ 



tan \ 7 



D 



(A4) 
(A5) 



where L2 is the radius of the cavity and L\ is the radius 
of the contact. For the small values of the ratio L1/L2 
this mode is also uniform in the most of the volume. In 
this case the three dimensional mode becomes 



g(r) oc 



Li _Li 
2r 



Li 



L? 



(A6) 
(A7) 



7 = ZL X DIL\ 

and its shape is similar to that which is shown in Fig. |^ 



APPENDIX B: THE OPTIMAL FLUCTUATION 
METHOD IN THREE DIMENSIONS. 

According to the optimal fluctuation method the prob- 
ability to, remain is given by the following system of 
equations^ 

p(t) oce 2 « , (Bl) 

t = irv I cosh(e)df, (B2) 
Js 



DV 2 9 = -iusmh(9) 



(B3) 



and the boundary conditions for the last equation were 
written after Eq._(pOh. This very problem was solved by 
Falko and EfetovEJTSec. VI and VII. In three dimensions 
one has 



r 



A/e J 



luJ 



3LiD 
A 



«1, 



7T 



p(t) oc e 



37r£> A*L X 
2A — TTT- 



«1 1„„2 tA 
2 = e ~2E [ °S — 



(B4) 

(B5) 
(B6) 

(B7) 



where 7 is given by Eq. (A7). This result is precisely 
equal to our one mode result Eq. (Q) in the long time 
limit. 
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